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I.  INTRODUCTION 

This  is  the  final  technical  report  for  Contract  No.  F30602-72 
C-0305 .  The  object  of  the  program  has  been  to  provide  theoretical 
studies  of  selected  topics  associated  with  light-beam  propagation 
through  a  turbulent  atmosphere.  The  efforts  were  in  direct  support 
of  the  programs  being  conducted  at  Rome  Air  Development  Center.  Among 
the  topics  considered  have  been  temporal  spectra  of  atmospherically 
induced  phase  difference  fluctuations,  an  analysis  of  the  angle  of 
arrival  apparatus  at  RADC,  and  a  study  of  the  Fourier  transform  imaqe 
restoration  procedure. 

It  is  the  object  of  the  present  report  to  summarize  these  studies 
with  emphasis  being  placed  on  certain  portions  of  the  image  restoration 
work  which  has  received  the  least  coverage  in  past  reports.  This  will 
now  be  done,  starting  with  the  angle-of-arri val  apparatus. 


II.  ANGLE-OF-ARRIVAl  APPARATUS 

The  angle-of-arri  val  apparatus  is  one  built  for  RADC  by  the 
Perkin-Elmer  Corporation  as  part  of  a  multi-purpose  apparatus.  Past 
descriptions  of  this  part  of  the  Perkin-Elmer  package  have  been  semi- 
quantitative.  The  object  in  this  description  was  to  obtain  a  more 
quantitative  understanding  of  the  device. 

Angle  of  arrival  measurements  are  in  general  useful  because  they 
give  angular  coordinates  of  objects.  One  generally  denotes  by  angle  of 
arrival  the  angles  defining  the  normal  to  the  wavefront  of  a  light 
wave  originating  from  a  point  object.  If  there  are  no  intervening 
index  fluctuations  then  the  wavefront  will  be  a  spherical  one  centered 
on  the  object  with  normals  indicating  the  direction  to  the  object.  If 
there  are  intervening  refractive  index  fluctuations,  then  the  arriving 
wavefront  may  be  different  from  spherical  or  may  be  spherical  with  a 
different  center  or  radius  of  curvature,  giving  a  false  measurement. 

The  device  senses  angle  of  arrival  fluctuations  by  examining  the 
position  of  the  focused  diffraction  limited  point  image.  The  trans¬ 
verse  displacement,  5,  is  measured  and  results  are  interpreted  in 
terms  of  angle  of  arrival,  a,  by  the  standard  relationship  a  =  5/d  where 
d  is  the  telescope  image  distance. 

All  angle  of  arrival  measurements  are  limited.  Even  in  the 
ideal  case  measurements  having  better  accuracy  than  the  classical 
angular  resolution  determined  by  the  input  aperture  of  the  optical 
system,  1.22  A/D  for  a  round  aperture,  are  difficult.  If  there  are  no 
refractive  index  fluctuations,  the  object  of  a  stationary  point  source 
is  a  stationary  diffract, on  pattern,  the  position  of  the  object  being 
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determined  by  the  center  of  the  pattern.  If  there  is  a  SDlierical  wavp 
front  with  center  modified  by  the  index  fluctuate is ,  then  the  dif¬ 
fraction  pattern  will  be  centered  at  a  different  spot,  the  amount  of 

thl  oK  If  Jh  proportio"a,  t?  the  apparent  angular  deviation  of 
diffSti™  i-f-^7e  1S  a  n?n-sPherical  wavefront,  then  the  standard 
riPnlnHin  1J1  t6d  Patcer"  is  broken  up,  the  amount  of  degradation 
depending  on  the  aperture  diameter  compared  with  the  mean  transverse 

spot  DOS^ion  anHP  fl  uctuat!°"s  •  The  device  determines  the  mean 
spot  position  and  averages  over  the  spot  shape. 

The  angle  of  arrival  apparatus  is  shown  schematically  in  Fig.  1. 
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Fig.  1.  Angle  of  arrival  apparatus. 


cXrjr-:  sopr?e  by  the  turt>"'*"t  atmosphere  is 

lil  ted  soot  th  hr“irn9  In  eSCuPe  a]d  1ma9ed  t0  a  de9™ded  diffract 
mi  ted  spot.  The  spot  is  then  chopped  with  the  scannir.q  reticle  whee' 

"rod?  ’s  onW  detrt0r-  Ihe  d?teCtPr  output  has  the  ?«™  of  a 

su?f?c  ent  ?f  st  so  9'  V  ■  ^  scanni"9  rate  of  the  “heel  ^ 

tl  n  oi  v  1  so  that  the  spot  is  chopped  many  times  before  the 

phere  changes  it  appreciably.  The  detector  output  is  then  filter 
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Fig.  2.  Typical  detector  output. 


to  pass  the  fundamental  frequency.  This  filtered  signal  is  compared 
for  phase  with  a  standard  signal  derived  from  the  opposite  side  of  the 
same  scanning  wheel.  The  phase  difference  is  normalized  to  zero  mean, 
squared  and  averaged.  The  mean  square  phase  is  then  related  in  the 
analysis  to  the  angle  of  arrival  variance. 

The  results  of  the  analysis  are  presented  in  terms  of  one  com¬ 
ponent,  <x,  of  the  normal  ~  to  the  plane  best  fitted  to  the  incoming 
wavefront.  The  results  are  based  on  the  assumptions  of  a  Kolmogorov 
refractive  index  spatial  spectrum  and  the  assumption  of  negligible 
amplitude  effects.  The  expression  for  the  mean  square  angle  of  arrival, 
<a£>,  which  the  device  measures,  is  then  expressed  in  the  form 


2>  =  <<j>  =  VD/2ro) 


-At  =  -A - A  F(N  ) 

|k|2  2k2 (D/2 )2  0 


where 


i 


Mr)  =  6.88(r/rJ5/3  = 


X  =  2*/k 
Nq  =  k6D/2d 


atmospheric  wave  structure  function 
neglecting  inner  and  outer  scale 
effects 

receiving  telescope  aperture  diameter 


=  wavelength  of  incoming  light 

=  Fresnel  number  based  of  the  size 
of  the  slits  in  the  scanning  wheel 

=  width  of  square  wave  slits  in  the 
scanning  wheel 

=  image  distance  of  receiving 
telescope. 
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Equation  (1)  is  very  similar  to  one  given  by  Hufnaqel  (Hufnagel ,  1963) 
The  difference  is  the  slowly  varying  correction  represented  by  the 
function  F(N0)  of  the  non-dimensional  parameter  NQ.  Equation  (1)  can 
be  derived  by  rion-dimensional izing  the  pertinent  integrals.  The 
function  F(N0)  is  obtained  from  the  evaluation  of  a  six-dimensional 
computer  evaluated  integral.  It  is  shown  in  Fig.  3.  For  the  RADC 


FRESNEL  NUMBER,  N0 


Fig.  3.  Correction  factor  F(No). 


telescope  and  apparatus,  5  =  2.8  mm,  D  =  40  cm  and  10.6  micron  radiation 
N0  =  46.  so  F(N0)  =  1.56. 


Equation  (1)  is  the  first  result.  Other  pertinent  points  covered 
previously  (Collins,  1972)  are 

1)  The  device  does  not  measure  what  has  been  commonly  defined 
as  angle  of  arrival,  but  does  measure  something  with  similar 
aperture  symmetry  properties. 

2)  For  very  large  and  very  small  apertures,  arrival  angle 
variance  is  independent  of  N0  =  k5 D/2d  (5  =  scanning  wheel 
slot  size,  d  =  image  distance). 

3)  Only  light  entering  the  aperture  in  a  ring  of  width  dX/6  at 
the  aperture  edge  is  effective  in  determining  the  arrival 
angle  variance. 

III.  PHASE  DIFFERENCE  TEMPORAL  SPECTRA 

The  second  area  of  investigation  was  phase  difference  temporal 
spectra.  The  work  on  temporal  spectra  is  of  interest  because  it 
provides  information  on  time  scales  of  phase  and  arrival  angle  useful 
to  systems  designers,  and  further  because  phase  difference  fluctuations 
are  intimately  I’elated  to  angle  of  arrival  fluctuations  for  small 
apertures  (Zi  nts  master ,  1971). 

More  generally,  temporal  spectral  measurements  provide  a  simple 
method  for  manipulating  experimental  data  to  give  a  check  on  our 
knowledge  of  the  light-beam-turbulence  interaction.  It  is  much  easier 
to  instrument  time-spectral  measurements  because  data  need  be  taken 
for  only  one  apparatus  setting.  The  corresnonding  spatial  correlations 
or  spatial  spectra  require  many  measurements,  each  with  the  apparatus 
adjusted  for  a  different  separation. 

Temporal  correlations  and  spectra  are  in  themselves  useful 
quantities  to  compare  with  theoretical  estimates  to  gain  information 
about  atmospheri c- turbulence  parameters.  In  this  context,  they  do 
not  suffer  from  aperture  limitations,  as  do  spatial  measurements, 
because  readily  available  time  lags  are  large  enough  to  represent 
spatial  separations  much  larger  than  those  conveniently  available. 

Thus  they  provide  a  convenient  optical  means  with  which  to  examine 
outer-scale  effects,  which  are  important  in  the  low  temporal  frequency 
region  of  the  phase  difference  spectrum.  Further  we  will  show  that 
the  high  temporal  frequency  dependence  of  the  phase  difference  spectrum 
is  easily  related  to  the  small  separation  dependence  of  the  phase 
structure  function,  thus  allowing  a  convenient  experimental  check  of 

Ds(p). 
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The  physical  situation  considered  is  shown  in  Fig.  4,  where  a 


y 


Fig.  4.  Diagram  of  the  phase  difference  spectra  situation. 


wind  of  constant  velocity  V-j  is  blowing.  The  spherically  diverging 
beam  is  centered  on  the  z  axis.  The  +x  direction  is  defined  so  that 
Vj  =  V  x  +  V,|  z,  where  V  is  the  magnitude  of  the  transverse  wind  and 
V„  is  the  component  of  the  wind  parallel  to  the  propagation  path.  V(| 
may  be  either  a  positive  or  a  negative  quantity.  The  phase  difference 
(<t>2-<t>l)  of  two  phase  points  separated  by  a  distance  p  in_an  x-y  plane 
is  recorded  as  a  function  of  time.  The  vector  distance  p  makes  an 
arbitrary  angle  e  with  the  +x  axis,  thus  enabling  us  to  determine  the 
phase  difference  spectrum  for  all  relative  orientations  between  the 
phase  points  and  the  transverse  component  of  the  wind  (Vx).  Not 
shown  is  that  the  optical  axis  specifying  the  direction  of  the  light 
beam  may  be  either  vertical  or  horizontal. 

Typical  results  for  a  horizontal  path  are  presented  in  Fig.  5 
where  we  see  the  power  spectral  density  normalized  to 

8.77k2LC3p5^3  versus  ft,  frequency  normalized  to  2irp/V 

where  X  =  2Tt/k  =  light  wavelength 

L  =  range 

=  refractive  index  structure  constant 
p  =  receiver  point  separation 

V  =  wind  velocity. 


normalized  frequency 


Fig.  5.  Normalized  phase  difference  spectrum  for  a  horizontal 
path  versus  normalized  frequency  and  orientation 
angle,  when  the  normalized  separation  is  0.1. 


A  "Von  Karmann"  index  spatial  spectrum, 


•  GO  =  .033 


( 

V 


A"1 1/6 

) 


was  chosen  so  as  to  include  effects  of  the  outer  scale,  L0.  The  more 
general  case  where  the  exponent  (-11/6)  was  replaced  by  an  arbitrary 
power,  v,  was  also  considered. 


Figure  5  shows  the  general  form  of  the  frequency  spectrum  to  be 
expected  for  constant  horizontal  wind  velocity,  various  angTes  of  tilt 
of  the  observation  points,  and  point  separation  0.1Lo.  Also  shown  are 
lines  giving  asymptotic  expressions  extracted  from  the  mathematical 
expressions. 


In  many  situations  the  wind  will  have  vertical  components 
necessitating  an  averaging  over  angle  and  washing  out  of  detailed 
angular  differences  shown  in  Fig.  5. 

In  the  case  of  slant  paths  cjj  is  proportional  to  height  to  the 
(-4/3)  power.  The  previous  spectrum  then  is  multiplied  by  an  additional 
factor  depending  on  the  function  h  of  the  maximum  and  minimum  heights, 
^max  anc^  Hmin, 


(2) 


h  =  (H  -  H  .  )/H  .  , 
max  mm  mm 


The  additional  factor  is  graphed  in  Fig.  6  as  a  function  of  h  for 
various  relative  frequencies.  The  case  considered  is  a  downward  slant 
path.  The  upward  slant  path  case  was  considered  also  and  gives 
different  results  (Collins,  1972). 

This  work  was  prematurely  terminated  to  allow  for  the  study  of 
atmospheric  image  restoration  to  be  considered  in  the  next  section. 


IV.  IMAGE  RESTORATION  WORK 

The  final  area  of  investigation  to  be  considered  is  a  theoretical 
examination  of  the  Fourier  transform  method  of  restoration  of 
atmospherically  degraded  images.  Ihe  incentive  for  this  work  came 
from  the  need  to  provide  information  for  a  proposed  series  of  experi¬ 
ments  on  the  usefulness  of  the  Fourier  transform  restoration  procedure 
for  restoring  satellite  images  using  an  airplane-born  reference 
source.  The  work  described  answered  the  question  whether  a  collimated 
beam  would  be  better  than  a  point  source  on  an  aircraft  flying  at  a 
much  lower  altitude  than  the  satellite.  The  formalism  is  also  useful 
in  predicting  the  isoplanatic  patch  size  for  this  restoration  scheme. 

The  Fourier  transform  restoration  procedure  is  well  known 
(Harris,  1966).  There,  one  images  simultaneously  a  turbulence  degraded 
object  and  a  turbulence  degraded  point  reference  source.  If  every 
point  in  the  object  and  the  point  reference  source  are  degraded 
exactly  the  same,  then  the  quotient  of  the  degraded  object  spatial 
spectrum  and  the  reference  point  spatial  spectrum  is  the  restored  image 
spatial  spectrum.  This  is  the  ideal  case.  The  investigation  reported 
herein  considers  the  situation  where  the  reference  source  is  different 
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Fig.  6.  The  additional  factor  for  downward  slant 
path  phase  difference  spectra. 


in  location  from  the  object  so  that  its  image  has  slightly  different 
degradation  from  that  of  the  object. 

The  situation  of  interest  is  shown  in  Fig.  7.  We  consider  two 
point  sources,  possibly  at  different  distances  from  an  imaging 
telescope.  The  light  from  the  two  sources  travels  through  a  turbulent 
atmosphere  and  is  imaged  on  two  image  planes  appropriate  to  their 
distances.  One  then  records  the  instantaneous  intensity  patterns  and 
uses  them  with  the  appropriate  transverse  size  scaling  in  the  Fourier 
transform  restoration  procedure. 
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Fig.  7.  Schematic  illustration  of  physical  situation. 


The  restoration  procedure  in  this  situation  is  expected  to  have 
limitations  The  situation  becomes  complicated  if  (1)  tne  object  po 
anT reference  point  are  not  at  the  same  transverse  location  so  that 
ra  s  f“  ?wo  partly  intersect  areas  of  different 
tions  and  (2)  the  object  and  point  reference  sources  are  at  different 
distances  from  the  reaver  so  that  refractive  Index  i f  uctua  ions  o  a 
given  size  experience  different  magnifications  by  the  time  they  have 

arrived  at  the  receiver. 

One  way  to  overcome  the  limitation  arising  when  the  object  is  at 
infinity  and  the  reference  source  is  at  a  finite  distance  from  the 
receiver  is  to  use  a  collimated  rather  than  a  Point  reference  source 
This  wolild  duplicate  the  point  source  at  infinite  distance  and  should 
work  if  the  degrading  turbulence  is  between  the*  reference  source  a 

the  receiver. 

The  object  of  these  calculations  is  to  assess  the 
displacement  of  the  reference  source  and  to  determine  the  feasibi  ^ 
of  using  the  plane-wave  source  as  a  possible  remedy  in  overcoming  the 
effects  of  longitudinal  distance  difference. 

A  model  is  chosen  which  allows  the  mathematics  to  be  greatly 
simplified.  The  model  is  written  in  terms  of  the  de^5tlon^  °^he 
turbulence-degraded  light  wavefront  from  the  unde9raded  wavefront 
These  deviations  are  described  mathematically  in  terms  of  a  ser  es  o 
polynomials  in  transverse  position,  r,  which  indexes  the  deviatm  at 
anv  aiven  aperture  position.  Thus  we  represent  the  phase,  ♦,  of  the 
right  arriving  at  potion  r  in  the  aperture  from  a  given  source  in 

the  form 

(3a)  <t>(r)  =  <t>0(r)  +  ^(r) 


where 


(3b)  -Mr)  =  l  a  F  (r) 

n 

and  where  ♦„  represents  the  phase  In  the  ideal  (non-turbulent)  case  and 
*i  (r)  represents  the  atmospherically  Induced  fluctuations.  The  poly 
JimUl  series  similar  to  that  used  by  Fried  (Fried.  1966)  are  given 
in  Eqs .  (4),  up  through  second  order. 


/3  _x_ 
a  a/2 


« 


(4c) 

(4c! ) 

(4e) 
(4f ) 


f_  =  —  ~y~ 

3  a  a/2 

1  —  —  I  1,2  t  J2-  !C\  \  I  t  -  tn\  2 


2  \  2  (x  +y  -  (a  /6))/(a/2)‘ 


F5  -  H\||  (»V)/(a/2)2 

f6  =  (3/a)(xy)/(a/2)2. 


These  polynomials  are  chosen  so  as  to  be  orthonormal  over  a  square 
aperture  of  side  "a".  It  is  assured  that  the  major  effect  of  the 

borSuH.ri^r56"^11  by  the  firSt  Six  ?olyno,,rials-  an  assumption 
by™r?ed  1  1  the  mean  square  sense^  by  our  calculations  and 

Drorpdurp  1°  thl  prob]em  has  two  parts  describing  restoration 

procedure  and  Jie  atmospheric  interaction,  both  of  which  can  be  stated 

the  ^ftitites  in  Eqs .  (3).  In  the  first  or  restoration 
part,  the  Fourier  transform  restoration  procedure  is  carried  out  usinq 

In^he1™^  thdel  f0e  the  ob:iect  and  reference  input  field  wavefronts. 

In  the  model  the  wavefronts  are  deterministic  and  are  represented  bv 

side^H  rttermS  1P  t?6  series,in  Eq.  (3b).  Indeed  for  the  case  con¬ 
sidered  later,  one  polynomial  (the  same  one)  is  used  for  both  the  obiect 

and  reference  beam  phases.  The  two  wavefronts  differ  in  the  size  of  the 
coefficients  an  and  an'  multiplying  the  polynomials;  i.e.,  they  have  the 
same  shaped  wavefront  deviation  but  one  is  more  bent  than  the  other 
[he  quality  °f  the  restoration  is  then  investigated  as  a  function  of  the 

t zl the+c°eIIicients.  This  allows  a  limit  to  be  set  on  the 
down  °f  the  coefflcient  difference  before  the  restoration  procedure  breaks 


In  the.  second,  or  atmospheric  part  of  the  study,  the  coefficient 
difference  is  further  related  to  the  atmospheric  parameters.  Specifically, 
the  root -mean -square  coefficient  difference  is  calculated  as  a  function  of 
object-to-reference  source  separation,  range,  wavelength,  turbulence 
i™Zqth’  and.receivmg  plane  aperture  size,  etc.  The  combination  of  the 
image  processing  and  atmospheric  calculations  for  the  coefficient  difference 
then  allows  one  to  calculate  a  measure  of  restoration  success  for  S 
set  of  atmospheric  parameters.  y 

The  restoration  and  atmospheric  portions  will  now  be  outlined 
and  results  given. 


In  the  restoration  portion  of  the  studies  the 
of  the  object  and  reference  fields  were  obtained. 


correlation  functions 
These  then  represent 
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the  associated  image  spatial  spectra  and  the  quotient  of  the  two 
spatial  spectra  represents  the  restored  image  spectrum  from  which  the 
restored  image  was  obtained.  The  restored  image  was  then  examined 
using  three  measures  of  image  quality:  the  restored  image  strehl  ratio, 
the  integral  scale  width,  and  the  mean-square  difference  from  the  ideal 
image.  The  whole  process  was  performed  for  a  series  of  increasing 
differences  between  the  object  and  reference  wavefront  polynomial 
coefficients. 

Typical  results  of  such  a  series  of  calculations  are  shown  in 
Fig.  8.  There  we  see  the  normalized  integral  scale  width  as  a  function 
of  coefficient  difference  for  a  sphprical  distortion  of  0.75  wave¬ 
lengths  maximum  deviation  of  the  object  wavefront.  We  see  that  the 
normalized  restored  image  width  is  unity  for  zero  difference  between 
the  wavefronts  and  initially  decreases  with  wavefront  difference. 

There  are  further  extremely  large  values,  both  positive  and  negative. 
The  large  values  come  from  dividing  by  very  small  numbers  when  the 
spatial  spectrum  quotient  was  obtained  in  the  restoration  procedure. 
This  process  gives  rise  to  unnaturally  large  values  in  the  restored 
image  spectrum.  When  the  large  values  cancel  there  is  reasonable 
restoration.  When  they  do  not  cancel  the  restoration  procedure  fails, 
as  indicated  by  the  very  large  or  very  small  integral  scale  width 
values. 

One  conclusion  drawn  from  these  calculations  is  that  the  restor¬ 
ation  procedure  is  reliably  valid  only  within  a  very  small  range  of 
differences  between  the  object  and  reference  wavefronts.  For  larger 
wavefront  differences,  either  further  image  processing  must  be  used  to 
prevent  the  unnaturally  larqe  values  of  the  restored  image  spectrum  or 
a  probabilistic  description  must  be  used. 

The  second,  or  atmospheric  part  of  the  study  was  then  aimed  at 
calculating  the  root-mean -square  wavefront  deviation  arising  as  a 
function  of  atmospheric  and  other  physical  parameters.  This  is  then 
to  be  compared  with  the  deviation  allowed  on  the  basis  of  the  restor¬ 
ation  procedure.  The  comp  ete  description  of  atmospheric  effects  is 
quite  lengthy,  so  only  the  pertinent  steps  will  be  documented. 

Referring  to  Fig.  7  for  the  physical  situation,  wn  start  by 
denoting  by  (7%r-| )  the  wavefront  arriving  at  7  in  the  aperture 
originating  from  an  object  at  7j ,  and  by  an  its  wavefront  expansion 
coefficient.  Then  using  the  orthonormality  of  the  functions  Fn(r)  in 
conjunction  with  Eq.  (4b)  we  have 

(5)  an  =  dr  Fn(r)$ (r)W(r) 
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Graph  of  normalized  restored  image  integral  scale  width  versus  wavefront  differed 
for  object  wavefront  deviation  of  three  quarters  of  a  wavelength. 


where  W(r) 


1  -a/2  <  x  <  a/2,  -a/2  <  y  <  a/2 

0  otherwise 


Similarly,  denote  the  wavefront  from  the  source  at  T|  +  p"  by 
4> ('P‘,T'i  ,+p“) ,  and  denote  its  wavefront  expansion  by  a^  given  by  an 
expression  similar  to  Eq.  (5)  but  with  r-j  replaced  by  ri  +  p.  If  the 
two  sources  are  superimposed,  then  p"  =  0  and  an  =  a^.  For  finite  values 
of  p",  an  and  ap  are  different  and  a  measure  of  their  difference  is 


(6) 


Dan(B>  *  <an"an>^  ' 


where  the  brackets  indicate  ensemble  average.  From  Fq.  (5)  we  have 

2 


(7a) 


D  (p)  =  < 
an 


j  U-,  (r  ,r-| )  -  <t>1(r,r1+p))Fn(r)W(r)dr 


(7b) 


J  dr  dr1  Fn (r )Fn (r * )W( r )W ( r ' )D^(r,r '  ,p ) 


where 


(7c)  D  (r,r',P)  =  <(^1(r,r1+p)-^1(r,r1))(<fri  (r‘  ,r1+p)-^(r*  ,r]))> 


phase 


The  atmospheric  effects  are  introduced  in  the  expression  for  the 


(8) 


♦  (r,r-,+p)  =  k  n  (r3^z3 


and  in  the  expression  for  the  index  structure  function 
(9)  W?,)  =  <(n(ri)-n(rj))2>  =  ^(^(r.+rj)) 


rrrj 


2/3 
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Using  Eqs.  (7),  (8),  and  (9),  the  expression  for  D^(r,r'+p)  simplifies 


(10)  D  ( r , r+p  )  =  S- 


L  rL 


0  "0 


dzl  dz2Cn  (^’(zl+z2  >  >  x 


r5”r4 


2/3+ 


r3-r6 


2/3  -  - 
“  r3"r4 


2/3 


Vr6 


2/3 


where  r3**r6  represent  points  along  four  rays  connecting  two 
aperture  plane  points  to  the  two  source  points. 

A  general  expression  for  D^r/r'p-)  can  be  derived  in  quite  a 
straightforward  fashion.  It  involves  writing  the  equation?  for  the 
four  1 ines^connecting  points  r  and  r1  in  the  aperture  piano  with  source 
points  at  q  and  q+p",  (points  rV •  -Tg  respectively  are  on  these  lines) 
and  then  substituting  in  terms  of  sum  and  difference  coordinates 

(lla)  6  =  z2-z1 

(llb)  a  =  (z2+zi  )/2 

(llc)  X  =  r-r' 

(lid)  n  =  (r+p  )/2. 


The  final  general  expression  is  included  here. 


(12b) 


[(Cx-((p1+nx+^x)(o"r5)/Ll^(^2"nx+K)(o+^)/L2))2 

+  (^y-  ( (r|y+2^y )  )/L-| )+( (ny-^y )  (a+^<5  )/l_2  )^+5^]^ 


(12c) 


r3~r6 


2/3  _ 


[(V((p2-nx'K)(a'^)/L2)'(-pTT1x+K)(a4S)/Ll))2 

1/0 

+  (  V  ( ("y^y  >  I0-?5  )/L2 )+( (ly-^y)  >/L,  )  )  V] 


(12d) 


)/l2)-((p2-i/^x)(o+1«  )/l2)  )2 

♦<y(  <y?y  (»4*  )/l2)+«vK  >(o+^  )/l2  )  )2+«2]1/3 


(12s) 


2/3  = 


[(5x+((-pi-nx-2ix)(a-^)/L1)-((-p1-T1x+^x)(o+l6)/Ll  ^  ^ 

+( v ( (vK)(a^  )/Li  )+((vK)(a+^  )/Li }  )2+62j1/3- 
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Ir  tq<;.  (12)  and  Lj  are  the  distances  of  the  reference  source  and 

co'o?d?n!tesenf1thly  ^  P?  a'ld  P]  are  ro5Pectivoly  the  transverse 
coordinates  of  the  reference  source  and  the  object:  p-=pw?.  It  has 

further  been  assumed  that  the  lateral  source  displacements  are  only 

Fn  t/7i  \  direttlon-  Wlth  the  variable  transformations  in  tqs.  (11) 
tq.  I/bj  becomes  ' 


(13) 


D  (p)  = 
an 1 


d5'dnFn(n+^)Fn(ri-^)W(r1+^)W(n-^)D0(c,n,p). 

lennthlhe  °f  EpS  ‘  °2)  and  (13)  is  strai qhtforwd  but 

lengthy.  It  has  been  performed  for  three  special  cases:  obiect 

at  9ltIiT'?tCehb0t^at  anfl?^e  distance’  corresponding  to  an' object 
®  J1.1-1 e  height:  and  collimated  reference  beam  comina  from  an  air¬ 

craft,  object  at  infinite  distance  and  reference  at  finite  distance 
t°  obiect  at  satellite  heioht  and  point  reference  on  an 

to  Sc?ni;aan^n5JeH  3"d  reference . b<>th  at  finite  distance,  corresoondinq 
to  using  a  poin,  reference  on  an  aircraft  to  restore  the  image  of  another 

aircraft  at  che  same  height.  The  expressions  have  also  been  evaluated 

tn  ft2eJ1tU!r°n..lIh?*re  the  refractive  index  fluctuations  are  confined 
sniirro  d°  ta  function"  layer  at  a  given  height,  and  where  the 

source  and  reference  are  at  the  same  and  different  heights.  The  integrals 
were  evaluated  numerically.  yrdlb 

Hiffo  T^pical  rGsults  are  shown  in  Fiq.  9.  There  we  see  the  RMS  wavelength 
difference  corresponding  to  coefficient  a4  graphed  as  a  function  of  p  for 
point  sources  at  20  km  and  150  miles  and  point  reference  at  20  km.  Further 
two  turbulence  conditions  were  considered.  The  two  bottom  curves  show 
the  effects  of  surface  level  turbulence 


?  -4/3 

<Cn<z)  ’  <*/*„>  ) 


for  the  two  different  objects  and  the  two  top  curves  show  the  effects 
of  a  high  altitude  turbulence  layer  at  17  km  altitude.  These  curves 
apply  specifically  for  coefficient  a4  which  imnlies  only  a  spherical 
deformation  of  both  object  and  reference  beams.  Similar  curves  for 
coefficients  a^  and  ag  with  hyperbolic  wavefront  centered  on  the  axis 
and  aperture  diagonals  are  shown  in  Figs.  10  and  11  respectively. 

The  strenqth  of  the  high  altitude  turbulent  layer  was  chosen  so 
as  to  provide  the  same  values  for  a  plane  wave  wave-structure  function 
at  the  ground  for  a  point  source  at  infinite  range  and  turbulent  layer, 
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as  due  to  ground  lovol  turbulence ,  and  point  sourer  at  infinity.  Iho 
calculations  wore  made  for  \  -  0.5p  and  2.1  x  10"lb  at  one  meter 
height. 

The  curves  exhibit  some  expected  behavior.  For  the  cases  where 
point  object  and  reference  are  both  at  20  km  height,  the  RMS  wavefront 
difference  goes  to  zero  as  p  becomes  small.  This  is  expected  since 
light  from  both  (point)  sources  ends  up  traversing  regions  of  the  same 
refractive  index  fluctuations.  For  the  case  where  one  source  is  at  a 
different  height  than  the  other,  there  is  a  residual  value  for  the  RMS 
wavefront  difference  even  when  p  =  0.  This  is  because  a  given  index 
fluctuation  magnifies  light  differently  from  sources  at  different  dis¬ 
tances,  and  there  is  complete  cancellation  only  for  very  low  level 
turbulent  fluctuations.  It  is  also  noted  that  the  effect  of  the  high- 
level  turbulent  layer  is  significantly  greater  than  that  due  to  the 
ground  level  turbulence.  That  is,  a  small  separation  of  rays  from  the 
two  sources  produces  a  much  greater  affect  in  interacting  with  high 
level  index  fluctuations  than  with  ground  level  fluctuations,  due  to  the 
wavefront  magnification  of  spherical  waves. 

The  curves  can  be  useful  in  deciding  whether  restoration  is 
possible  in  a  given  situation.  For  example,  if  a  value  of  0.01 
wavelengths  difference  can  be  tolerated  for  the  wavefronts  then  a  very 
close  spacing  between  reference  and  image  must  be  obtained  or  there 
must  be  weaker  turbulence.  If  a  value  of  0.10  is  acceptable,  then, 
with  the  value  of  chosen,  a  source  separation  of  one  aperture  size 
is  allowed. 

One  other  point  is  indicated.  That  concerns  the  relative  ease 
of  restoration  for  objects  at  the  same  height  as  the  reference, 
compared  with  objects  much  higher  than  the  reference.  Figure  3  indi¬ 
cates  that  for  object  and  reference  separated  by  more  than  one  aperture 
size,  the  ease  of  restoration  is  almost  the  same,  the  values  of  the 
RMS  wavefront  difference  lying  within  a  factor  of  two  of  each  other  for 
p/a  >_  1 . 

Further  work  remains  to  be  done  in  this  area  especially  on 
the  restoration  portion  of  the  effort.  That  involves  investigating 
the  use  of  a  smoothing  function  in  the  restored  image  spectrum  to 
remove  the  unnatural  singularities  coming  from  dividing  by  small 
numbers.  That  will  provide  more  physically  meaningful  restored  images 
and  better  numbers  on  the  restored  image  quality  indices  with  which  to 
compare  the  mean  square  wavefront  deviations. 

V.  GENERAL 

in  addition  to  the  work  outlined,  the  bibliography  lists  one 
paper  that  has  been  published  Reinhardt,  1972)  and  three  that  are  in  the 
preparation  stages;  (Davidson  1974),  (Collins,  1974a),  (Collins,  1974b). 
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VI.  SUMMARY 


This  contract  has  considered  several  problems,  all  aimed  at 
providing  support  to  work  being  carried  out  in  connection  with  RADC 
programs.  The  problems  are  phase  difference  temporal  spectra 
predictions,  analysis  of  the  angle  of  arrival  measurement  apparatus 
used  at  RADC,  and  theoretical  study  of  the  Fourier  transform  method 
of  restoration  of  atmospherically  degraded  optical  images.  In  each 
case  theoretical  predictions  were  completed  and  where  called  for 
recommendations  made  and  conclusions  drawn. 
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